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Abstract 

According to a fluid dynamic expansion of the fireball we investigate how the inhomogeneity 
of the system influences the chiral phase transition of QCD. We compare the averaged values 
of the order parameter in equilibrium with that of a homogeneous system. If the temperature 
is averaged over a certain region of the fireball the corresponding correlation length does not 
diverge in an expansion with a critical point. 



1. Introduction 

The search for a possible critical point of QCD is one of the driving motivations for the 
low energy runs at RHIC/BNL and the upcoming facilities FAIR at GSI/Darmstadt and NICA 
in Dubna. The assumed existence of a critical point is based on two observations: at vanish- 
ing baryochemical potential highly sophisticated lattice QCD calculations revealed an analytic 
crossover and at high baryochemical potentials various approaches to effective models in 
mean-field and beyond find a first order phase transition, which then necessarily terminates in a 
critical point Q S B m ■ 

At both, the first order phase transition and the critical point, an enhancement of fluctuations 
is predicted to be seen in heavy-ion collisions. For the first order phase transition these are based 
on the growth of fluctuations due to the nonequilibrium effect of supercooling leading to nucle- 
ation and spinodal decomposition flHSIstl- Such fluctuations should be visible in studies of 
single event particle distributions. The critical point, however, is supposed to imprint its signals 
in event-by-event fluctuations of conserved quantities, such as net-charge or net-baryon number 
distributions ifloifTlll . In thermodynamic systems the correlation length diverges. Systems cre- 
ated in a heavy-ion collision, however, are finite in size, highly dynamic and a temperature and 
baryochemical potential can at most be defined locally. The finite size of the system naturally 
limits the growth of the correlation length. The finite time, which the dynamic systems spends 
at a critical point, limits the growth even further because relaxation times also become infinite 
at the critical point. This phenomenon is called critical slowing down. Assuming a phenomeno- 
logical time evolution of the correlation length with parameters from the 3d Ising universality 
class it was found that the correlation length does not grow beyond 2 - 3 fm [12J. In all of 
these approaches the background medium is assumed to be homogeneous in its thermodynamic 
quantities, like the temperature. 

Recently, the model of nonequilibrium chiral fluid dynamics (N^FD) for the explicit prop- 
agation of fluctuations coupled to a dynamic expansion of a heavy-ion collision was developed 
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Il3ll showing that nonequilibrium effects can have an important influence on the evolution of the 
fluctuations of the order parameter of chiral symmetry (1141 IX 511 . In this model the expansion of 
the fireball created in a heavy-ion collision is described by ideal fluid dynamics, where the fluid 
consists of quark and antiquarks with the local thermodynamic properties of the quark-meson 
model |2|| . 

In this work we want to investigate the influence of the inhomogeneity of the fluid dynamic 
background on the sigma field and the correlation length. For this purpose we restrict ourself to 
a local equilibrium treatment of the sigma field, instead of an explicit propagation but work with 
concepts that have already been studied in NjfFD. We refer to these references 

M El for 

more information. 



2. Equilibrium chiral fluid dynamics 

For a fixed baryochemical potential the strength of the phase transition can be tuned via 
the quark-meson coupling g. At hb — and for a realistic coupling of g — 3.3 the transition is a 
crossover, for g = 3.63 it is a critical point with T c a 140 MeV and a first order phase transition 
with T c a 123 MeV is found for g = 5.5. 

The thermodynamic potential of the quark-meson model in mean-field approximation is 

a(T) = -^\nZ=U(cr,7f) + a q g(T), (1) 

with the chiral potential U (cr, n) and the quark contribution to fermionic one-loop level. The 
equilibrium values for the chiral fields, <x eq and 7? eq , are obtained from minimizing the thermody- 
namic potential (Q} with respect to cr and n. While cr eq becomes finite at lower temperatures due 
to dynamic symmetry breaking, n vanishes for all temperatures. We do not continue to explicitly 
mention the pion field in the notation. The minima of the thermodynamic potential O define the 
thermodynamically stable states of the matter. In this case the pressure and the energy density of 
the system are 

p(T) = -QCOU^ , e(T) = T \^p-j ~ P(T) ■ (2) 

The fermions are described fluid dynamically by the relativistic equation of energy-momentum 
conservation 

dyj^ = , (3) 

with the energy-momentum tensor of an ideal fluid. Starting from initial conditions from UrQMD 
lfl6Tl the fluid has an inhomogeneous temperature distribution, which is determined by the evo- 
lution of the energy density. Therefore, the equilibrium value of the sigma field is defined only 
locally in each cell of the fluid dynamic grid. Speaking of the temperature of the system thus im- 
plies some form of averaging over parts of the system. As a consequence, all quantities deduced 
from thermodynamic expressions are subject to the same averaging. 

A prominent feature in a scenario with a first order phase transition is the discontinuity of 
the sigma field at the transition temperature, which will be studied in section [3] In a critical 
point scenario we investigate the correlation of the fluctuations based on the definition via the 
curvature of the thermodynamic potential in section |4] 
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3. Discontinuity in a first order phase transition scenario 

In figure [1] we present the temperature and the sigma field averaged over an inner volume 
of about 30 fm 3 . As the system expands it dilutes and the energy density and the temperature 
become smaller until the coexistence region of the first order phase transition is reached. Here, 
the local temperature stays constant. Since almost all of the volume has an energy density in the 
coexistence region, also the averaged temperature stays almost constant. The averaged sigma 
field increases during the transition. If we assumed a homogeneous temperature of the system 
with the same evolution as the averaged temperature in figure Q] and assumed the corresponding 
equilibrium value of the sigma field for the entire system, we recover the discontinuity in the 
evolution of the sigma field. This discontinuity is washed out when averaging over an inhomo- 
geneous system. 
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Figure I : Scenario with a first order phase transition. Time evolution of the averaged temperature (solid line) and the 
averaged sigma field (dashed line) in a fluid dynamically expanding system, where the local sigma field is in equilibrium 
with the fluid. Assuming that this temperature evolution was the global temperature of a homogeneous system, we show 
the global sigma field (dotted line) for comparison. 



4. Correlation length in a critical point scenario 

In the right plot of figure |2]the same curves are plotted as in figure[T]but here for a scenario 
with a critical point. Again we find that in the transition region the evolution of the sigma field 
is washed out in an inhomogeneous system after averaging compared to an assumed homoge- 
neous temperature evolution. In the left plot we show the influence of the inhomogeneity on the 
correlation length, which is locally evaluated from 

i - *?m (4) 

We observe that due to averaging the correlation length in an inhomogeneous system does not di- 
verge at the averaged critical temperature of the system. It remains finite and grows up to 1 .5 fm. 
When we assume again that the temperature evolution was a homogeneous temperature of the 
system we can calculate a single correlation length, which would correspond to the system's 
correlation length. It is shown in figure |2]for comparison. 

In finite systems one speaks about a critical region around the critical point, where the corre- 
lation length is still significantly enhanced but does not diverge. For an inhomogeneous system 
we find the same situation when averaging over parts of the system. 
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Figure 2: Scenario with a critical point. The curves in the right plot are the same as in figure^but for a critical point. 
The left plot shows the averaged correlation length (solid line) compared to a correlation length corresponding to a 
homogeneous system temperature (dashed line). 



5. Summary 

We showed that the inhomogeneity of a system created in a heavy-ion collision smears out 
the evolution of the order parameter of chiral symmetry, the sigma field, in the transition region. 
The effect on the correlation length is more dramatic, as its divergence at the critical temperature 
in a homogeneous system is reduced to an increase up to 1 .5 fm. In the investigated scenario, 
however, fluctuations of the sigma field are not directly studied. It needs to be said that the growth 
of the correlation length is a very dynamic process and should be studied in a fully dynamic setup, 
what is currently worked on within N^FD. 

This work was supported by the Hessian Excellence Inititive LOEWE through the Helmholtz 
International Center for FAIR. We are grateful to the Center for Scientific Computing (CSC) at 
Frankfurt for providing the computing resources. 
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